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EQUATIONS FOR THE RAPID MACHINE COMPUTATION 
OF EQUILIBRIUM COMPOSITION OF AIR AND DERIVATIVES 
FOR FLOW- FIELD CALCULATIONS 

By G. Louis Smith, Wayne D. Erickson, 
and Mary R. Eastwood 
Langley Research Center 

SUMMARY 

Equations are presented for the rapid computation of equilibrium air composition 
for use in flow-field calculations by electronic computer. The temperature-density 
domain with temperatures from 1500° K to 15 000° K and densities from 10^ to 10“ 7 
times standard atmospheric density is covered by four equilibrium air models which 
describe the composition in terms of mole fraction to an accuracy of better than 0.001. 
Equations are also presented for the computation of partial derivatives of mole fractions 
with respect to temperature, pressure, and elemental composition. 

INTRODUCTION 

The rapid computation of the chemical equilibrium composition of air is of great 
importance in flow-field problems in which a very large number of equilibrium composi- 
tion calculations must be carried out. The importance of extreme rapidity in this com- 
putation is enhanced because, in a Runge-Kutta numerical integration over a single 
interval, it is necessary to compute the equilibrium composition at each of several steps 
within the Runge-Kutta procedure. The computation of chemical equilibrium composition 
may be accomplished by a number of methods. The methods of Brinkley (ref. 1), Huff, 
Gordon, and Morrell (ref. 2), and White, Johnson, and Dantzig (ref. 3) appear to be the 
most widely used general methods. A comparative study by Zeleznik and Gordon (ref. 4) 
showed that none of these methods had any significant computational advantage over the 
other two. In reference 5, Erickson, Kemper, and Allison show that by reducing the sys- 
tem of equations for a particular chemical system to one equation in one unknown, the 
unit computational speed may be increased by two orders of magnitude over these three 
general methods. This gain in computational time is at the expense of tedious algebra 
in reducing the system to a single equation. 



In flow problems with equilibrium chemistry, the partial derivatives of mole (or 
mass) fractions with respect to pressure and temperature are needed. If diffusion is 
being considered, so that the elemental composition is variable also, the partial deriva- 
tives of mole (or mass) fractions with respect to elemental composition are also needed. 

The purpose of this report is to present the application of the method of refer- 
ence 5 to air over a wide range of temperature and pressure. In this paper, the tempera- 
ture range from 1500° K to 15 000° K with the density range from 10^ to 10" ^ times 
standard atmospheric density is divided into four regions. Within each region a model 
of air and the applicable composition equations are set up. Following the method of ref- 
erence 5, the system of equations for each air model is reduced to one equation in one 
unknown mole fraction. This equation may be rapidly solved numerically. The 
remaining concentrations may then be computed from explicit relations. In addition to 
basic equations for computing composition, this paper also presents these partial deriva- 
tives of composition. 


SYMBOLS 

A function of defined by equation (41) 

a m function of elemental composition defined by equations (31), where sub- 

script m is any integer 

B function of x^ defined by equation (42) 

bjj ratio of number of species k atoms to number of nitrogen atoms, where 

subscript k is 2 for oxygen and 3 for argon 

D determinant for equations (37) and (38) 

d n coefficients in quartic equations, defined by equations (15), (62), and (77) 

for models A, C, and D, respectively 


F 1 

function of 

x 2 

defined by equation (A2) 

F 8 

function of 

x 7 

defined by equation (A22) 

F 11 

function of 

x 7 

defined by equation (A31) 
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partial pressure equilibrium constant for reaction j; subscripts for j 
are given in table II 

boundary between models 

mean molecular weight of mixture 

molecular weight of species i 

integer 

pressure, atmospheres 
Universal gas constant 
temperature, degrees Kelvin 

temperature behind normal shock wave, degrees Kelvin 
mole fraction of species i; subscripts for i are given in table II 
functions of xj defined by equations (33) 
functions of xj defined by equations (35) 

functions of elemental composition defined by equations (16), where sub- 
script m is any integer 

change of enthalpy across reaction j at standard atmospheric pressure 
ratio of density to standard atmospheric density 

mole fraction equilibrium constant for reaction j ; subscripts for j are 
given in table II 

ANALYSIS 


The basic approach of dividing the density- temperature domain into four regions, 
each with its appropriate air model, is first discussed. Next, for each air model, the 
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equations for equilibrium composition are set up, the method of solution is briefly indi- 
cated, and the final equations required for the determination of equilibrium composition 
are presented. The limitations of each model are then indicated. Finally, the equa- 
tions for the partial derivatives of composition with respect to pressure, temperature, 
and composition are presented, with an indication of the method of derivation. 

The composition of air is described to a 0.1-percent accuracy with respect to the 
tables of Hilsenrath and Klein (ref. 6) which were "taken as the standard for composition 
for equilibrium air. The present report does not include gas imperfections, which are 
considered in reference 6. In figure 1, the air density for which the mole fraction of 
any given species equals 0.001 is shown as a function of temperature. Thus, for an 
accuracy of 0.001, a species is present in significant amounts on one side of one of its 
lines but not on the other side. The symbol for the species appears on the side of the 
line for which the species is present in amounts greater than a mole fraction of 0.001. 
After an inspection of figure 1, it was decided to represent the composition of air in dif- 
ferent density-temperature regions with four appropriate models as listed in table I. 

The region of validity of a model is the region within which every species present in a 
mole fraction of 0.001 or greater is included in the model and is indicated in figure 1 for 
each of the four models by heavy lines. It is seen that together the four models cover 



Figure 1.- Lines for mole fraction equal to 0.001 and regions of validity of air models. 
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the density- temperature domain covered by reference 6 (i.e., 1500° K ^ T 2 15 000° K; 

-7 = log 2 2.2 or 2.4) except under conditions of very high density, where species such 

as NC >2 and N 2 + become significant. These regions not covered by the four models are 
not important to the aerodynamicist, inasmuch as they occur only in the stagnation 
region of a blunt body at velocities exceeding 20 000 ft/sec (6.1 km/sec) at altitudes 
less than 40 000 ft (12.2 km). 

The species are numbered as indicated in table II. Also listed in table II are the 
nine independent reactions and the corresponding equilibrium constants in terms of mole 
fractions. 

A discussion of the different air models follows. 

Model A 

Model A is the low temperature model and allows for dissociation and NO forma- 
tion. The region of validity is limited on the high temperature boundary by the line for 
electron mole fraction equal to 0.001 and on the high density boundary by the line for 
NO 2 mole fraction equal to 0.001. In this high density region, gas imperfections are 
becoming significant, which further limits the accuracy of this model on the high density 
boundary. Below 1500° K, model A is still valid, but the only significant departure from 
cold air composition is due to the formation of NO. No computer computations were 
made for temperatures below 1500° K, but no difficulties are anticipated. 

Model A consists of the six species listed in table I: The six equations for 
model A involving the mole fractions of these species are 


T 1 = X 4/ X l 2 

(1) 

t 2 = x s/ x 2 2 

(2) 

r 3 = x 6y / Xl x 2 

(3) 

x 3 /( x l + 2 x 4 + Xg) = b 3 

( 4 ) 

( x 2 + 2x 5 + x 6)/( x l + 2x 4 + x 6) = b 2 

( 5 ) 

X 1 + x 2 + x 3 + x 4 + x 5 + x g = ! 

(6) 


Equations (1), (2), and (3) are equilibrium conditions, equations (4) and (5) are mass 
constraints, and equation (6) follows from definition of mole fraction. Equations (1) to 
(4) may be used to express Xg, x^, Xg, and Xg in terms of x^ and Xg as 
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X 4 = T lXl 2 

(V) 

x 5 * r 2 x 2 2 

(8) 

x 6 = r 3 x l x 2 

(9) 

x 3 = b 3 (! + 2 ti x i + t 3 x 2 )xi 

(10) 


Equations (7) to (10) are used to eliminate x g through x g in equations (5) and (6). 
This procedure yields two equations in Xj and x 2 ; thus, equations (5) and (6) become, 
respectively, 

x 2 + 2r 2 x 2 2 + r 3 x l x 2 = b 2( x l + 2r l x l 2 + r 3 x l x 2) (“> 

and 


x 1 +x 2 + t^Xj + 2 Ti*i 2 + t 3 x i x 2) + T l x l 2 + t 2 x 2 2 + r 3 x l x 2 = 1 (12) 

Further manipulation yields 


and, finally, 


where 


2b 2 ~ y 0 x 2 ~ y 4 x 2 2 

b 2 + n x 2 


l 

n-f 


d x 9 n = 0 
n l 


d 0= b 2 2 (^ + 4 


di = -bo 4y ft + 


r 2 j_ri 
T l 


d 2 = ( r 0 2 ~ 4 b 2^) " 


h 2 y 3 + y l y 2 




d 3 * 2 r 4 v 0 - b 2 5 y, - 'Tl 


2 T 2 

4 = y \ “ y l T{ 


d A = - 2 


(13) 


(14) 


(15) 
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( 16 ) 


r 0 = 1 + 2 ( b 2 + b 3 ) 

r l = t 3 ( r 0 " b 2) 
y 2 = 1 + b 2 + b 3 + r 3 (l - b 2 ) 

y 3 = t 2 (2 + b 2 + 2bg) + T 3 ( b 2 + b 3 ) 

^4 = 2t 2(^0 “ b 2) 

Equation (14) is in the form suggested t>y reference 5 and can he solved quite rapidly t>y 
using Newton’s iteration scheme. Once x 2 is found, can be calculated from equa- 
tion (13), and the remaining species are computed by equations (7) to (10). 

Equation (14) has four roots. Only real roots lying between 0 and 1 need be con- 
sidered. As pointed out in reference 5, for only one of these values of x 2 will the mole 
fractions of the remaining five species lie between 0 and 1. 

Model B 

Model B is the most complex of the air models. Covering a wide range of tempera- 
tures, model B allows for dissociation, ionization, and formation of NO. The anticipated 
region of validity of this model was limited by the line for 0 2 mole fraction equal to 0.001 
on the low temperature boundary, by the line for A + mole fraction equal to 0.001 on the 
high temperature boundary, and by the NO + and N 2 + mole fractions equal to 0.001 on the 
high density boundary. Unfortunately the full extent of the anticipated region of validity 
was not realized due to numerical difficulties, which are subsequently discussed. As a 
consequence, the practical region of validity as limited by these numerical difficulties 
was found by trial and error to be bounded by the line for electron mole fraction equal 
to 0.001 on the low temperature boundary. This limitation creates no problem in prac- 
tice, inasmuch as model A can be used up to this line. 

Model B consists of the eight species listed in table I. The eight starting equations 
for this model are 


x 4 /xi 2 = 

(17) 

x 6/ x l x 2 = t 3 

(18) 

x l/ x 7 x 8 = t 4 

(19) 

x 2 / X 7 X 9 = t 5 

(20) 

x 3/( x l + x 8 + 2x 4 + x 6) =b 3 

(21) 

( x 2 + x 9 + x e)/( x l + x 8 + 2x 4 + x e) = b 2 

(22) 

x? = x 8 + x 9 

(23) 

x l + x 2 + x 3 + x 4 + x 6 + jc? + x 8 + Xg = 1 

(24) 
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Equations (17) to (20) are equilibrium conditions, equations (21) and (22) are mass con- 
straints, and equation (23) expresses charge neutrality. Equations (17), (18), and (21) 
may be solved for x 4 , Xg, and Xg, respectively, and these unknowns are eliminated 
from the system of equations by the following substitutions: 


X 4 = TjXj- 


x 6 “ t 3 x 1 x 2 


x 3 = b 3(xi + x 8 + 2x 4 + x 6 ) (27; 

Equation (22) is used to eliminate Xg; thus, 

Xg = b 2 (xj + x 8 + 2r lXl 2 + TgXjXg) - x 2 - TgXjXg (28) 

Equation (23) is now combined with equation (28) to give the following expression for Xg 

x g = (1 + b 2 )" 1 [x 7 + x 2 + t 3 x iX2 - b 2 ( X j + 2r lXl 2 + TgXjXg] (29) 

Next x 7 is solved for from equation (24) with equations (25) to (29) used to eliminate 
Xg, x 4 , Xg, x 8 , and Xg. After some manipulation the result is obtained that 

2 

x 7 = a 0 ” a l x l " a l x 2 " a 2 T 3 x l x 2 “ a 2 r l x l (30) 


a o = ( 2 + r^ 


a l = a-of 1 + TTbi 


a 2 = a ol 1 + ITbi 


For convenience equation (30) is written as 

x 7 = Oi - 0^2 


a \ ~ a 0 “ a l x l " a 2 r l x l 


“2 = a l + a 2 T 3 x l 




( 34 ) 


Equation (32) is used in equation (29), whence Xg may be written as 

x 8 = ^1 + % x 2 


where 


h = (* + b 2) _1 [ a o " ( a l + b 2) x i “ ( a 2 + 2b 2) T i x i^ 1 

h= { 1+ b 2)" 1 [( 1 " a l) + i 1 ~ a 2 " b 2) T 3 X l] J 

> 

(35) 

Equation (23) gives, by using equations (32) and (34), 



x 9 = (“i “ %) " ("2 + ^) x 2 


(36) 

Thus far, all the species, in particular x^, Xg, and Xg, can be expressed as very 
simple functions of Xj and Xg, which remain undetermined, and equations (19) and (20) 
have yet to be used. Equations (32), (34), and (36) are used in equations (20) and (21). 
After some manipulation equations (20) and (21) become 

-( 2a l“2 + Tj) x 2 + X 1 - a l 


(37) 

(“A - Vl) x 2 - “2^2 X 2 2 = ^ x 1 ' Vl 


(38) 


Equations (37) and (38) may be considered to be two equations in two unknowns 
therefore, solving by Kramer's rule gives 

x = ^~ 

2 D/cfcj 


where 


A = + a 1 2 /3 2 - ^■ x l( a 2 + %) 

(41) 

B M ( 2 “l“2 + £t)(“A - 7J x l) + (“A - Vl)(“l 2 - ^ x l) 

(42) 

(°i“2 + ^2 + “2% 

(43) 

Equations (39) and (40) give X£ and Xg 2 as functions of Xp which is as yet undeter- 
mined. The equation for x^ is obtained by squaring the right-hand side of equation (39) 


X2 and 

(39) 

(40) 
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■I nil II I ■■■■him II anil 


and setting it equal to the right-hand side of equation (40), whence, after simplification, 


/ D\_ A S 

^j B - “2 A 


Equation (44) has the form F(xjJ = 0 and may be solved numerically for x^, after 
which the remaining Xj may be computed. In principle, equation (44) may be expanded 
to give an eleventh- order polynomial; however, in practice, the tedious labor required to 
carry out this expansion makes it prohibitive to do so, so that equation (44) is used 
directly. 

It was found that for low temperatures where electron concentrations are small, 
equation (30) is unsuitable for computing x ? after x 1 has been determined, inasmuch 
as a small number is obtained by subtracting large numbers with a resulting loss of 
accuracy. The same problem appears in equations (34) and (36). This problem is cir- 
cumvented by rewriting equations (19) and (20) as 


X 7 X 8 = ^ X 1 


X 7 X 9 = T^ X 2 

Adding these two equations and using equation (23) gives 


X7 = F 4 Xl + ? T 2 

Now Xg and x g can be computed from 


X S = r 4 x? 

x 9 = f§ (47 

y t 5 x 7 

Once equation (44) has been solved for x^, x 2 follows immediately from equa- 
tion (39). Next x 4 , x 6 , x ? , x g , and Xg are found from equations (25), (26), (45), 
(46), and (47), and finally x g is found from equation (27). 


Model C 


Model C applies at high temperatures and consists of atoms, ions, and electrons. 
It is limited on one boundary by the presence of N 2 and on the other boundary by the 
presence of N ++ . 
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Model C contains the seven species listed in table I. The equations used for this 
model are 


x l/ x 7 x 8 = t 4 

(48) 

x 2 / X 7 X 9 = Tg 

(49) 

X 3/ X 7 X 10 " t 6 

(50) 

( x 2 + x 9)/( x l + x 8) =b 2 

(51) 

( x 3 + x 10)/( x l +x 8) =b 3 

(52) 

x 7 = x 8 + x 9 + x 1Q 

(53) 

+ x 2 + x 3 + x 7 + x 8 + x 9 + x i 0 = 1 

(54) 


Equations (48), (49), and (50) are equilibrium conditions, equations (51) and (52) are mass 
constraints, and equation (53) expresses charge neutrality. Equations (48), (49), and 


and Xg to give 


X 1 = r 4 x 7 x 8 

(55) 

x 2 = t 5 x 7 x 9 

(56) 

x 3 = r 6 x 7 x 10 

(57) 


By using equations (55), (56), and (57), equations (51) and (52), respectively, give 


= b 2 1 + TjX, X 8 


1 + t 4 x 7 


x 1n = b 


1 + T 


4*7 


10 ~ "3 i" 


o 


+ t 6 X 7 

Equations (58) and (59) are used in equation (53), whence 


(58) 

(59) 


1 + t 4 x 7 

1 + b 9 - + b 

2 1 + r 5 x 7 


1 + t 4 x ? 
3 1 + r 6*7 


( 60 ) 
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Thus, only x^ remains to be determined; equation (54) is the only equation left unused. 
Equations (55) to (60) are used to express all xj in equation (54) in terms of x^. The 
result, after simplification, is 

4 

£ d nX7 n = 0 (61) 

n=0 


where 

d Q = -(1 + b 2 + b 3 ) 

dj = 2(1 + b 2 + b 3 ) - (b 2 + b 3 )r 4 - (b 3 + 1 )t 5 - (l + b 2 )r 6 
d 2 = (l + 2b 2 + 2b 3 )r 4 + (2 + b 2 + 2b 3 )r 5 + (2 + 2b2 + b 3 )Tg 

- ( t 5 t 6 + b 2 T 6 T 4 + b 3 T 4 T 5) 

d 3 = (i + b 2 + 2b 3 )r 4 r 5 + (2 + b 2 + b 3 )r 5 T 6 + (l + 2b 2 + b 3 )r 6 r 4 
d 4 = (1 + b 2 + b 3 )r 4 r 5 T 6 


Equation (61) can be solved by Newtonian iteration to give x^. Then, Xg is computed 
from equation (60), Xg and x^q from equations (58) and (59), and x^, x 2 , and x 3 
from equations (55), (56), and (57). 


Model D 

Model D applies at extremely high temperature where the air consists of singly 
and doubly ionized atoms and electrons. It is limited on one boundary by the presence 
of N; the other boundaries were not considered. 

Model D and its equations are analogous to model C and its equations. The atoms 
and singly ionized atoms of model C correspond directly to the singly and doubly ionized 
atoms of model D. The equations of model D are 


y vn - t 7 

(63) 

Wia ■ t 8 

(64) 

*10/ X 7 X 13 ■ r 9 

(65) 
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( x 9 + x 12)/( x 8 + x ll) =b 2 

(66) 

( x 10 + x 13)/( x 8 + x ll) " b 3 

(67) 

x 7 = x 8 + Xg + X 1Q + 2(x n + x 12 + x 13 ) 

(68) 

x? + x 8 + ^ + x io + x n + x 12 + x 13 = 1 

(69) 

These equations are required by equilibrium conditions, mass constraints, and charge 
neutrality. 

Solving equations (63), (64), and (65) for Xg, Xq, and x 1Q 

results in 

x 8 = r 7 x 7 x n 

(70) 

x 9 = t 8 x 7 x 12 

(71) 

x 10 = t 9 x 7 x 13 

(72) 


Equations (70) to (72) are used in equations (66) and (67) to give 


x 12 

x 13 


- b 1 * , 

D 2 1 + TgX, X xi 

1 + T„Xn 

= b 3 1 + T x Xl1 
1 + t 9 x 7 


(73) 

(74) 


Equations (70) to (74) are used in equation (68) to give the following relation between x 7 
and 


X 11 “ 


*7 

A + t 7 x 7\ 

2 + r 7 x, + bJ rr; p-j(2 + TgX,) + b< 


(l + Tr i Xr n 
A + t 9 x 7/ 


(75) 


f2 + r, 


9 * 7 ) 


Equations (70) to (75) are used to express all xj in equation (69) in terms of Xij. The 
result simplifies to 

4 

£ d n x 7 n = 0 (76) 

n=0 
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where 


^0 - " 2(1 + b 2 + ^ 3 ) 

dj = 3(1 + b 2 + b 3 ) - t ? ( 1 + 2b 2 + 2b 3 ) - r g (2 + b 2 + 2b 3 ) - r g (2 + 2b 2 + b 3 ) 

d 2 = r 7( 2 + 3b 2 + 3b 3) + r 8( 3 + 2b 2 + 3b 3 ) + T g( 3 + 3b 2 + 2b 3 ) 

- t 7 t 8 (1 + b 2 + 2b 3 ) - t 8 t 9 (2 + b 2 + b 3 ) - T g T 7 (l + 2b 2 + b 3 ) 

d 3 = t 7 t 8( 2 + 2 ^2 + 3 ^3) + t 8 t 9( 3 + 2 ^2 + 2b 3 ) + r 9 r 7 ( 2 + 3 ^2 + 2b 3 ) 

- t 7 t 8 t 9 (1 + b 2 + b 3 ) 

d 4 = 2 t 7 t 8 t 9 (1 + b 2 + b 3 ) 


> ( 77 ) 


Equation (76) is solved for x 7 by Newtonian iteration, after which is computed 

from equation (75), xj 3 an< l x j 3 f rom equations (73) and (74), and x g , x g , and x^q 
from equations (70), (71), and (72). 


Partial Derivatives 

The partial derivatives of composition with respect to pressure and temperature 
are now considered. These derivatives enter into a flow problem through the equation 
of state 

y M 

which is expressed in derivative form as 

i VP = i Vp + I VT + MV(^) 

The term v (j^) is given by 

v ( s )=- pZ M ‘ Vx ‘ 

i 

For an equilibrium gas 

Vx ‘ = (If)t Vp+ (w) p vt 

/9xi\ / 9x-j\ 

In computing and (-— j , the pressure and temperature dependencies are 

implicit in Ty The approach taken in this analysis is to compute the values of 9xj/9-7j, 
after which the desired partial derivatives are given by 
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where the summations are carried out over only those j values for which r j appears 
in the gas model. All Tj are defined so that 


Thus, 


v»t; p v \ st / p rt 2 j 


The equation in one unknown Xj of the form 

F(Xi,Tj) . 0 


is differentiated implicitly to give 


ax T 


8F/8Tj 


dri 9 F/3xj 


The remaining 9x^8 Tj are derived by differentiating the explicit expressions for each 
Xi . The derivations are straightforward, and the resulting expressions are listed in 
appendix A in the order required for computation. 


The partial derivatives of mole fractions with respect to elemental composition 
Sx i/ 9b k are needed for the solution of problems involving diffusion, for example, in 
boundary layers. Again, the derivations are straightforward, and the resulting expres- 
sions are listed in appendix B in the order required for computation. 


Minor Species 

A simple computation of the concentration of a minor species can be made in which 
the major species are known from the solution of one of the air models. The method is 
demonstrated by considering the following example. 

Assume that the electron concentration in the region of model A is desired. Elec- 
trons are contributed primarily by NO but are also contributed by N and O. The 
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equations applicable to the problem for the indicated chemical reactions are as follows: 


N + + e = N 

x l/ x 7 x 8 ’ t 4 

(78) 

0 + + e = O 

X 2/ X 7 X 9 - t 5 

(79) 

NO + + e = NO 

X 6/Vl4 - T 10 

(80) 


x 7 = x 8 + *9 + x 14 

(81) 


In equations (78) to (80), xj, x 2 , and x 6 are known from the major species computa- 
tion and Xrj, Xg, x g , and x J4 are unknown minor species. Solving equations (78), 
(79), and (80) for x g , Xg, and x 14 results in: 




X 8 t.x, 


4 7 


= t 5 x 7 


(82) 


x 14 = 


x 6 


r 10 x 7^ 


These expressions are substituted in equation (81) to give the following equation for x^: 


x 7 = 



x 2 ^6 
t 5 t 10 


(83) 


The electron concentration is thus quickly computed from equation (83). If desired, the 
mole fractions of the minor species NO + , N + , and 0 + can then be computed from 
equations (82). 


RESULTS AND DISCUSSION 


The equations for the four gas models were programed for machine computation, 
and a number of computations were made to test the accuracy of the gas models in com- 
puting composition. No numerical difficulties were encountered with models A, C, or D, 
but problems were encountered with model B at low temperatures, which effectively 
limit its region of validity. At very low temperature, x 4 becomes quite small as Tj, 
t 4 , and Tg become quite large compared with unity, and in equation (44) the terms 

IL b and OjA 2 are insensitive to Xj and are equal to each other to several decimal 
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places. Practically, equation (44) reduces to an identity and is satisfied quite accu- 
rately by a fairly wide range of x^; this creates a computational problem. Another 
problem encountered in using model B at low temperatures is that an extraneous root of 
equation (44) is obtained; however, Xg computed from this value for x^ exceeds the 
bound imposed by the mass balance of the problem so that ^ can be tested, this root 
of Xj rejected, and the search for a satisfactory root continued. None of these prob- 
lems arise with model B on the high temperature side of the line for electron mole frac- 
tion equal to 0.001. (See fig. 1.) Inasmuch as model A is applicable out to this line, all 
the numerical difficulties associated with model B at low temperatures can be avoided by 
using model A up to the line for electron mole fraction equal to 0.001. The results of a 
series of test computations for air composition are listed in table III; also given for com- 
parison are the compositions from reference 6. This comparison indicates that the 
accuracy of the composition of air as computed from the equations in the present report 
should be adequate for most purposes. 

The boundaries of the gas models, shown in figure 1 as a density-temperature plot, 
were mapped onto a velocity- altitude plot (fig. 2) by using figure 6 of reference 7. Thus, 
the temperature and density data in figure 2 correspond to equilibrium conditions imme- 
diately behind a normal shock wave. Compositions are listed in reference 6 for tempera- 
tures up to 15 000° K. The boundaries below the 15 000° K isotherm were extended out 
to 50 000 ft/sec (15 km/sec) by using the compositions listed in reference 8. The low 
density boundary ^log^Q = -7^ is above an altitude of 400 000 ft (122 km). It is not 



Figure 2.- Velocity-altitude plot of regions of validity of air models. (Data based on conditions 
immediately behind normal shock.) 
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likely that equilibrium flow conditions will be of concern above this altitude. The high 
density region where NO + and Ng* are significant maps onto a velocity -altitude region 
corresponding to orbital velocities below a 40 000-ft (12 km) altitude. These flight con- 
ditions are not currently important to the aerodynamicist. 

The pressure-temperature boundaries of the four models are shown in figure 3. 
For purposes of machine computation, quadratic expressions were fitted to the bound- 
aries. The boundaries established by these expressions and indicated in figure 3 by 
bold lines are: 

I: T = 1500° K 
II: T g 15 000° K 

mt log 10 p = -4.6 - 0.035(l5 - - 0.00556(15 - ^ f 

IV: lo gl()P = 3 

V: l = 1 + - 8) - 0.16f~~ - 8) 2 

V1000 / V1000 ) 

for log^Q p = Z, use model A 

for logjQ p < l , use model B, C, or D 

V I: l *= 0.5 + 0.446^ - 12 ) - 0.0541 (jT. - lif 
for log^Q p > l, use model B 
for log 10 P = i, use model C or D 

V”: 1 - - 2 ' 65 + 0 ' 51 (i^0 ' 15 ) 

for log^Q P = l, use model C 
for logjQ p < l, use model D 

The species listed for model A are adequate for T < 1500° K, but the equations have not 
been applied in this region to check for numerical difficulties which may arise. Simi- 
larly, no check has been made of the models at temperatures above 15 000° K or at pres- 
sures below boundary m. Thus, boundaries I, n, and III are included to avoid using the 
models in untested regions and are not tested bounds on the validity of the models. How- 
ever, boundary m corresponds to altitudes in excess of 400 000 ft (122 km), and for tem- 
peratures below 1500° K the composition of air is essentially that of cold air. Bound- 
ary V was placed so as to avoid computational problems with model B by using model A 
as much as possible, and boundary VI was placed so as to use model C as much as 
possible, because model C is considerably simpler than model B. 
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Figure 3 .- Pressure-temperature boundaries of air models. 


RESUME 

The computation of air composition has been dealt with by dividing the temperature 
range between 1500° K and 15 000° K and the density range between 10^ and 10“^ times 
standard atmospheric density into four regions, within each of which a relatively simple 
air model may be used for the equilibrium chemical composition. In each air model, 
the mass constraints and equilibrium conditions, stated as a function of temperature and 
pressure, are reduced to a single equation in one unknown mole fraction. This equation 
can be solved numerically, after which the composition of the gas can be rapidly com- 
puted from closed-form expressions. The accuracy of the compositions computed by the 
present method was better than 0.001 in terms of mole fraction. 

Equations are also given for partial derivatives of mole fraction with respect to 
temperature, pressure, and elemental composition. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., April 27, 1967, 

124-07-01-29-23. 
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APPENDIX A 

8x i 

EQUATIONS FOR COMPUTING 

The equations are presented in the order required for computation for each model. 
For model A, 

4 


3*2 n 4 8T i 2 

8r< 4 


The egressions for 8d n /8rj for model A are given in table IV. The equations used in 
the computation of 8x^/8 tj for air model A are as follows: 

. *>2 - rft - W m 

b 2 + V2 

SFi y 0 + 2 y 4 x 2 + rfr 

8 x 2 b 2 + n x 2 


8 Fj 

8x 2 


8x 2 

8t 1 


9 Fj 

8x 2 _ 

2 ( y 0 " b 2) x 2 2 

8x 2 

8t 2 

t>2 + yi ^2 

tH 

CD 

8x 2 _ 

(r 0 - b 2 ) Xl x 2 

8x 2 

8t 3 ■ 

b 2 + n*2 


93^ 9Xi o 

87J- ' 2l l x l S7J- + x i 


— = 2t,x, 

9r 2 * 8 t 2 


8x 4 8Xj 

wj - 2t 1 x 1 
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8Xr 8X 9 

= 2t 9 x 9 -z-Z 
3 Tj ^ ^ 3 tj 

8 Xr 9x„ 

9^ - 2t 2 x 2 9^ 

3x 5 3 x 2 

3 = 2t 2 X 2 9^ 



(a) 

(b) > (A6) 


(c) 


+ x 

3ri 3y)ri 2 1 QTjy 

!^ =T te x +x ^£\ 

&T2 3^9t2 2 1 9t2 j 


9 xg 3 y 9 Tg 2 1 873^ X 1 X 2 


(a) 

(b) > (A7) 

(c) 

J 



3 1 0Ti 


+ 2 


54 

STi 



(A8) 


The equations for air model B are 

9A _ 2„ e n- m a2 + 2b 2 - 2 

3 Tl " 2 X 1 “2^1 “l^ 1 + b2 X 1 


2 1 - a 2 - b 2 


9^ " Wl*! + “1 T7 


b 2 


x l“ 


8A 

875 


= 0 


"\ 


- 2a^/3j a 2 Xj2 

(a) 

2 L + 1 ' a 2"M 

(b) 

T 4 1 \ 2 1 + b2 / 


(c) 


(d) 


> (A9) 


3^ " “ 2a 2 x l 2Q! 2( Q! l^l ■ ^ x l) " ( 2a l a 2 + 7g)(Vl + a l ■ f —b^j x l 2 
~ 2 («i% " «2^l) a 2 x l 2 “l 

wj 

(Equations continued on next page) 
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9B 
9r 3 


“ 2a iVi( a ih - tj) x i + *i(“i 2 - x i)(“i 1 i* 2 ^ 2 ' a 2^i 
" ( 2 “l“2 + “1% ‘ “2% + 7^ 


9B _ 1 

9 t 5 


T^l^l^i " 74 x l 


a 2 + 2b 2 1 


Th) _ .. 2 a L e , - 

9r x " X 1 a 2^2 + “2 j + b2 


’( D /“2). 2a a s x +(2a a. , 

9 t 3 ” 20 ! l a 2 ^ 1 X l + + T5 j 1 + b2 X 1 


' 1 “ a 2 ” b 2 \ / \ 

“1 — iYb, ' a 2^lJ«2 x l - (“^2 - Vl] a 2 x l 


!m)_ 0 

9 t 4 

8 ( D / Q ^) % 

9rc ^ 2 


9a?-, 

_ = - a l - 2a 2 r lXl 


^-ar 
9 X1 " a 2 T 3 


8 ^1 _ a l + b 2 2( a 2 + 2b 2) T l x 

9xj 1 + b2 1 + 1>2 1 


9^ ^ 1 ~ a 2 - b 2 
9xj 1 + b£ 


(f) 

(g) 

(h) 

(i) 

> (A9) 

(j) 

(k) 

( l ) 

J 

(a) " 

(b) 

> (A10) 

(c) 

(d) 
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8A 9«1 da 2 9 0i „ 9a 1 9 ^9 , 

^ = ^^l + ^^^ + a l a 2^ + 2ffll ^^2^ 1 2 ^-^^ 2 ) 


H(^2 + ^2\ 

t 4 \ 9x i 9 *iy 


“1 ST - u) + (Sf % * “1 SI - § h - “2 1^)(“1 2 - ^ *1 


\/ a «1 


9x l 


♦(Vs- vOKsj-^- 


■ 2 (^- “2 + “1 5^2 + ( 2 “l“2 + 7^ - (®J H * °1 

80!;, 8j3,\ 8a, 

- 8^1 " <*2 iS^p “ (“1^2 " “2^l) 


(a) 


> (All) 


(b) 


(c) 


da 2 

- — = a^x, 
9t 3 2 1 




9«2 

9tT 


8F 

dTz 


= 0 


_ K 0 /^) B , p 8B 8Qi 2 A 2 _ 2Q! A ^A 

B + 0> 8t, 8t 4 ^ “ 

^ ] 3 


9T 3 


Bt a 


9F _ 8 ( P / Q! 2) b . D 8B _ 0a 2 a 2 _ 2 A 8A 

8 Xl - 8 Xl B + a^ A 


(a) 

(for j * 3) (b) 

(c) 

(d) 


> (A12) 


8x x 9F / 9r j 

8t^ _ " 8F/8 Xl 


(A13) 


23 



APPENDIX A 


9A 9A 8x l 
9 x 2 ^ 9Tj + 9 Xl St 7 

8t 3 * D /^ 



(A14) 


9 x 7 9 X 2 ^^1 2 

9r l ” 8x l 8t 1 " ^ 8t 1 8x l 8t 1 2 2 1 


(a) 


9x 7 9 q!j 9xj 9x 2 9 9xj 

9 Tg 9xj 9Tg ^2 9x3 9 xj 9T3 X 2 a 2 X l X 2 

9 x 7 9 q!j 9 x^ 9 x 2 9 Og 9 xj 

9r^ ” 9 x-£ 9t^ °^2 9r^ 9x^ 9r^ X 2 


(b) 

> (A15) 

(c) 


9x 7 9q!j 9Xj 9x 2 9q^ 9x^ 

9 ^“ 9 ^ 9 ^ _a! 2 9 ^“ 9 ^ 9^ X2 


(d) 


= !£l!fl + fl 8x 2 , 8 % 8x 1 a 2 + 2b 2 2 

9tj 9xj 9t^ 2 9t^ 9x^ 9tj 1 + bg 1 


ffft 


0x 8_ 9 ^1 8x 1 fl 8x 2 , 8 ^2 8x l , 1 - a 2 ' b 2 , y 
9t 3 9x 1 9 t 3 + 9t 0 9Xi 9r 0 1 + b„ 12 


= — - — + /3« 


8t 5 8 Xl 9 t 5 


8t 3 

9 xj 

8t 3 

8x 2 


9 x^ 

8t 4 + 

9 Xj 


8x 2 


9 Xj 

8t 5 + 

9X 1 

8t 5 


(a) 


(b) 


> (A 16 ) 


(c) 


(d) 


^9 = ^ 7_^8 

9 t. 9 r. 9 t. 

] 3 3 


(A 17 ) 


24 



APPENDIX A 


9x 4 

9 -T 

2T 1 X 1 

9x l 

9T 1 


2tiX. 

9x l 

9t 3 

1 1 

9t 3 

! X 4 = 

2r 1 x 1 

9x l 

9t 4 

1 1 

9t 4 

5i_ 

9t 5 

2t 1 x 1 

9X 1 

9t 5 


(a) 


(b) 


> (A 18) 


(c) 


(d) 


!fe = r +x 

9t 1 2 \ 9t 1 2 1 9r lj 

^- T f^l x + x !^\ + xx 
8t 3 " T 3^r 3 X 2 + X 1 drj + X 1 X 2 

!^_ T ( 8 _h x + x 

9t 4 3^8 T4 x 2 + X 1 9 T4 j 


(a) 


(b) 


> (A19) 


(c) 


!^_ T fe x +x 

8 T5 - ^ 9T5 X 2 + X 1 er 5 J 


(d) 




8Xg 

8T4 


+ 2 


8r. 


+ 



(A20) 


For model C, 



9t- 4 

£ nd n x 7 n_1 


n=0 


(A21) 


The egressions for 
used in the computation of 


9 dn/ 8 Tj for model C are given in table V. 
9xj/8rj for air model C are as follows: 


The equations 
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1 + t 4 x 7 1 + 
Fo = 1 + b 9 ±-1 + b„ 

° 2 1 + TgX 7 3 1 + T, 


T 4*7 


6*7 


(A22) 



f b 2 

9Xrj 4) 


+ r 5 x 7 

9F 8 1 


b 2 

9t 4 " X 7' 

i 1 

+ t 5 x 7 

1 

00 

Pn 

CD 

1 

+ 7-4X7 


1 + r 6 x 7 


[1 + r 4 x x 


T 5 b 2 T 6 b 3 

+ 


li 1 + T B x 7) ( 1 + r 6 x 7)' 




9F 8 


t 1 + t 5 x 7) 


2 b 2 x 7 


1 + r 4 x 7 


( 


1 + r 0 x. 


6 A 7 )' 


b 3 x 7 


(a) 

(b) 

(c) 


> (A23) 


(d) 




j./, 8F 8\ 8x 7 fa.f£a 

8-j Fj, \ 8 dXrjJdTj Fg STj 


(A24) 


3 x 


1 


874 



+ x 8 



+ X 8 X 7 






t 4 



x 7 



9x 9 _ b 2 T 4 x 8 ~ r 5 x 9 
9x 7 1 + t 5 x 7 

8x 9 _ ^9 _^7 ^ 8x 8 b 2 x 7 x 8 

3t 4 9x 7 9 t 4 + x 8 ,9 t 4 + 1 + T 4*7 

8x 9 _ 9 x 9 9 x 7 ^ x 9 9 x 8 x 7 x 9 

9t 5 9x 7 9t 5 x 8 ar 5 1 + 

^ 9 = ^ 9^7 ^ 9^8 

9t 6 9x 7 9t 6 + x 8 9t 6 


(a) 


Cb) > (A25) 


(c) 

J 


(a) 


(b) 


> (A26) 


(c) 


(d) 
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9x 


10 


9x 


10 


9r c 


9x<j 


8x 10 

= V 

4 X 8 

“ r 6 x : 

9x 7 


1 + 

rex? 

. 8x 10 

9x ? 

* + 

x 10 

8x 8 

dXrj 

8t 4 

x 8 

9t 4 

. 8x io 

9x 7 

x 10 

9x 8 

9 x 7 

8t 5 + 

x 8 

8t 5 

8x io 

!fl + 

x 10 

ff8_ 

‘ 8x 7 

8t 6 

x 8 

8t 6~ 

/ 9X, 

7 


8x 10 N 


: x io 

1 

+ x 7 

8t 4 j 


X 7 X 10 
+ r 6 x 7 


9x, 


'9 x ? 


9x 


9^" T6 l^; X 10 + X 7'^ 


10 ' 

5 ) 


9x 3 /9x 7 9x 10 \ 

9t 0 = T 6yarg X 10 + x 7 ~dTQ~J + X 7 X 10 


+ x ^ 

9t 4 5 \9t 4 7 9 9t 4 J 


9x<i 


' 9Xfl 


9x„ 


= Tr 


■ Xr, + 


9T5 ’ 5\9t 5 + 9t 5 


+ x„x, 


7 A 9 


_2_ r rfa x +x ^7 

9t r 5l 9r ft X 7 + X 9 Br, 


9Xr» /ft 

r 6 “\ or 6 

For model D, 


6 1 


9x^ 

9r4 


n=0 ^ 


4 

I 

n=0 


J ndnx ? 11 " 1 


(a) 


(b) 

> (A27) 

(c) 


(d) 


(a) 

(b) V (A28) 

(c) 


(a) 

(b) > (A29) 

(c) 


(A30) 


27 



Illll 


Him ii ii i in mini 


APPENDIX A 


The expressions for 9 d n / 9T j for model D are given in table VI. The equations 
used in the computation of 3x^9 Tj for air model D are 


b 2 b 3 


F H = ! + (! + t 7X7 )(1 + b 2 + b 3 ) + (1 + r 7 x 7 ) rv + ^ 


T x x 7 1 + t 9 x 7 j 


(A31) 


3Fii ( b« bq ) 

— — = r 7 1 + b 9 + bo + - — + ! — 

•\ * 3 1 + TgX 7 1 + T 9 X 7 1 


3 X 7 


" (1 + T 7 x 7)j 


b 2 T 8 + b 3 T 9 


(! + TqXtj) 2 (1 + TgXrjY 


3F 


11 


3rn 


= Xr 


1 + b 2 1 + 


l_\ + b 3 1 + — i- 

1 + t 8 x 7/ 3 V 1 + r 9 x 7/ 


9F, , 9: 


11 dx 7 


dXrj 


3F 


11 


9 To 


1 + TnX„ 3F, , 3x 7 
-b n Xr, “ 1 


2 7 


(! + r 8 x 7 ) 


2 9 Xrj 9 Tg 


3F 


11 


9r n 


1 + TrjXfj 9F 1:l 3x 7 


(a) 


(b) 


(c) 


(d) 


> (A3 2) 


3x 


11 1 


9F ll\ 9x 7 X 11 9F 11 


3Tj F n 


1 - x 


11 


9x 7 / 9T i F 11 9r i 


(A3 3) 


9 Xg 9 x 7 


9x 


11 


9^ ~ r 7 8^ X 11 + t 7 x 7 + X 7 X 11 


0X 8 _ 9x 7 t 9x 11 

8^~ 7 ^8 11 


r 8 


3x, 


9x 8_ t !Hl x +tx 

9t^ _ T 7 9t q X 11 + T 7 X 7 9 Tq 


(a) 


(b)> (A34) 


(c) 
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b 
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8t 7 ' 
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X 11 
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+ 
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9t 8 

3xiy 

9t 8 

X 11 
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!!i_ r +x !fi2\ 

a r 7 " 8^0r 7 12 + x 7 9r 7 j 


ax g _ / ax 7 3 x 12 \ 

9^ “ T8 \9^ Xl2 + X7 "9^7 + X?Xl2 


ax 9 /ax 7 ax 12 \ 

^=% x l2 + x 7^j 


(a) 


(b)> (A3 8) 


(c) 
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EQUATIONS FOR COMPUTING 


8b k 


The equations are presented in the order required for computation for each model. 
The following equations are used in the computation of 3Xy/3b k for air model A: 


8d 0 _ o d ° 

3b 2 b 2 


9d C 
3b c 


= 0 


8d l di , 

_A = _A- b 9 (l + 8 t,\ 

9b 2 b 2 2 \ 1) 


9di 

1E i = -b 2 (l + 8r 1 -2T 3 ) 


8d 2 r 

— = 4r 1 |2(b 2 r 2 - y oTl ) + 2b 2 (r 2 - 2^) + - y 3 - b 2 (r 2 + r 3 j 

- t 3?2 ~ Vli 1 + t 3) 


9d 2 

9b 3 

9d 3 

9b 2 


9dc 


= 4r l(^0 “ 4b 2 T l) " b 2( 2r 2 + t 3) “ 2t 3^2 “ y l 

= - 8T l( b 2 T 2 - n> T l) " 4 t 1^o( t 2 ' 2r l) - T 2 y l “ b 2 T 2 T 3 
" r 3^3 “ ^l( 2T 2 + T 3) 


— = -SrjTgbg - 2b 2 r 2 r 3 - 2T 3 y 3 - ri (2r 2 + r 3 ) 


0d 4 

— = Sr^bgTg - y 0 r 1 )( t 2 - 2 tj) - ^T 2 T 2 > y 1 
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8x2 

9b k 


I 


8d, 


n 


9b k 



4 

^ ndnXg 11 " 1 
n=0 


8 F i/axg is computed from equation (A2). 


a *1 aF l °*2 , „ 1 - x 2 - t 2 x 2 2 .. 
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1 + r 3 x 2 
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9x-» 


9F 1 9x 2 2 x 2 + 2 t 2 x 2" 
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8b 3 8x 2 9b 3 b 2 + n x 2 b 2 + n x 2 


8x 4 8x x 

8^ = 2 Vl 


8Xp; 9 x 9 
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0b k 2 2 9b k 


5« = t + x ^ 

9b k 3 l 9b k 2 1 9b k/ 


9x 3 

9b 2 

8x 0 


= b 


f'dx, dX A dXn) 
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9b 3 


9b 2 9b 2y 


/8x, 8x4 8 x r \ 
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The equations for air model B are 
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9b k \ 9b k ^ 1 9 b k / 2 \ 5jtt > k \ 8b k 2 1 9b k 
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dF/dx^ is obtained from equation (A 12 (d)). 
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The equations for air model C are 
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9do 

ibj - 2t 4 + 2t 5 + t 6 - t 4 t 5 


= r A r l 


'4 r 5 + t 5 t 6 + 2 t 6 t 4 


8d, 

8b 2 

8d3 

8b 3 4 5 5 


— - 2t 4 t 5 + t 5 t 6 + t 6 t 4 

u 

8d 4 

®bk 


= t 4 t 5 t 6 




n=0 


9d n n 

^ 7 


4 

I 

n=0 


J nd n x 7 n_1 


8F g 1 + t 4 x ? 

9b 2 1 + 

8Fg 1 + 

8b 3 1 + TqXtj 

9Fg/8x 7 is computed from equation (A23(a)). 



APPENDIX B 


Wi_ 8F 8\ ax 7 x 8 8F 8 

8b k F 8\ * 8 8x 7/ 8b k F 8 8b k 

8x 9 b 2 T 4 x 8 ~ r 5 x 9 8x 7 t 1 + T 4*7 L +b &X 8 \ 

8b 2 ~ I + T5X7 9b 2 1 + t 5 x 7\ 8 2 8b 2 / 


8x 9 b 2 T 4 x 8 ~ t 5 x 9 8x 7 , 1 + T 4*7 b f^8 

8b 3 “ 1 + 1-5X7 8bg 1 + T5X7 2 9b 3 

8x 10 _ b 3 T 4 x 8 ~ t 6 x 10 8x 7 t 1 * t 4 x 7 b &X 8 
8b 2 1 + tqXtj 9 b 2 1 + T 6 xrj 8 9 b 2 


8x 10 b 3 T 4 x 8 " t 6 x 10 8x 7 , 1 + t 4 x 7 ( h 
8 bq 


l + r 6 X7 8b 3 1 + t 6 x 7 


3 8b; + X 8 


8x l _ T te x +x 


8x 9 /8x q 8 x 7 \ 


8 X 3 

9b k 



8 x,\ 


The equations for air model D are 


8dQ 

8b 2 

8£o 

8b 3 

8dj 

~^2 


-2 

-2 


3 - 2t ? - Tg 


2 r n 



2 t„ 


2t 8 ' t 9 


8dn 

— = 3t 7 + 2t 8 + 3t 9 - t 7 t 8 


t 8 t 9 " 2 t 9 t 7 
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APPENDIX B 


a d 2 

^3 


= 3t 7 + 3 t 8 + 2 r g - 2t 7 t 8 - TgTg 


2t 9 t 7 


ado 

— i = 2r 7 r 8 + 2 t 8 t 9 + 3r g T 7 - r 7 r 8 r 9 


8<1q _ 

— = 3 r 7 r 8 + 2 TgTg + 2rgT 7 - 

9d 4 « 

= 2 tuTqTq 

abk 789 


Y n 

a „ Z/ ab k 7 

9x 7 n=0 k 


ab k 4 


^ nd n x 7 n 


n =0 


8Fn 

ab 2 

1 + r 7 x 7 

= 1 + r 8 x,( 2 + Vt) 




9F U 

= 1 + T ’ X ’(2*r 9X7 ) 




ab 3 

1 + TgX 7 V 9 1 ) 




9F ll/ 8x 7 

is computed from equation (A32(a)). 



8x ll 

1 L 9F ll\ 9x 7 X 11 8F H 




9b k 

" F ll\ ^7 ) 9b k F 11 9b k 




8x 12 

b 2 T 7 X ll ' t 8 X 12 9x 7 , 1 + t 7 X 7 / 

< 

bn 

9x n\ 

9b 2 

1 + TgX 7 9b 2 1 + t gx 7 ' 

v 11 

2 

ab 2 / 

9x 12 

b 2 T 7 x ll " t 8 x 12 8x 7 ( 1 + t 7 x 7 

bo 

axn 


9b 3 

1 + r 8 x 7 9b 3 1 + t 8 x 7 

2 

9b 3 


9x 13 

b 3 r 7 x ll " t 9 x 13 ^7 , 1 + r 7 X 7 

bo 

9x n 


ab 2 

1 + t q x 7 9b 2 1 + t q x 7 

u 3 

ab 2 




APPENDIX B 


9x 13 _ b 3 T 7 x ll - t 9 x 13 ^7 , 1 + T 7 X 7 / 8x ll 

Obg 1 + TgXiy 9bg 1 + TgXr j l 3 0bg 

!!§ = T /!5 x +x 

3b k 7 \®k 11 7 8b k / 

ex 9- r ^ x + x ^ 

^ V* 12 7 ®|J 
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TABLE I.- SPECIES IN THE AIR MODELS 
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TABLE H.- NUMBERING OF SPECIES AND REACTIONS 


Species 


Reactions 


Mole fraction 
equilibrium constant, 

T i 


i 


1 

2N = N 2 

T 1 = x 4/ x l 2 

1 

N 

2 

CM 

o 

II 

o 

CM 

t 2 = x 5/ x 2 2 

2 

O 

3 

N + O = NO 

t 3 " x 6/ x l x 2 

3 

A 

4 

N' 1 ' + e = N 

t 4 = x l/ x 7 x 8 

4 

N 2 

5 

0 + + e = O 

r 5 = X 2/ X 7 X 9 

5 

°2 

6 

A + + e = A 

t 6 " X 3/ X 7 X 10 

6 

NO 

7 

N ++ + e = N + 

t 7 = X 8/ X 7 X 11 

7 

e 

8 

0 ++ + e = 0 + 

r 8 = X 9/ X 7 X 12 

8 

N + 

9 

A ++ + e = A + 

t 9 = x lo/ x 7 x 13 

9 

o + 




10 

A + 




11 





12 

0 ++ 




13 

A ++ 




14 

NO + 
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TABLE m.- COMPOSITION COMPUTED BY PRESENT METHOD AND COMPOSITION FROM REFERENCE 6 


[First value of paired numbers, present method; second value, ref. 6] 


(a) Model A 


T, °K 

lo M^ 

N 

o 

A 

n 2 

02 

NO 

2000 

2 

0.2668 X 10" 10 

0.1026 X 10" 4 

0.9343 X 10" 2 

0.7770 

0.2054 

0.8233 x lO -2 



.2668 X 10" 10 

.1025 X lO' 4 

.9343 X lO’ 2 

.7767 

.2048 

.8219 x lO’ 2 


0 

.2916 X 10-9 

.1109 X lO -3 

.9343 X lO" 2 

.7769 

.2053 

.8274 x lO -2 



.2916 X 10-9 

.1109 X lO -3 

.9340 X lO -2 

.7767 

.2052 

.8270 x lO' 2 


-2 

.2918 X 10-8 

.1108 X lO' 2 

.9339 X lO" 2 

.7766 

.2047 

.8259 x lO' 2 



.2917 X 10-8 

.1108 X lO" 2 

.9335 X lO -2 

.7763 

.2047 

.8257 x lO’ 2 


-4 

.2903 X 10" 7 

.1090 X 10" 1 

.9292 X lO -2 

.7728 

.1989 

.8122 x lO' 2 



.2903 X 10 -7 

.1090 X 10 _1 

.9289 X lO" 2 

.7725 

.1988 

.8119 X lO’ 2 


-6 

.2785 X 10' 6 

.9270 x 10" 1 

.8911 X lO -2 

.7415 

.1500 

.6910 X lO -2 



.2784 X 10' 6 

.9269 X 10" 1 

.8906 X lO -2 

.7411 

.1500 

.6908 X lO -2 

4000 

1 

0.1210 X 10" 3 

0.4399 X 10" 1 

0.9137 + lO -2 

0.7166 

0.1357 

0.9447 x 10 _1 



.1210 X lO- 3 

.4394 X 10' 1 

.9135 X lO -2 

.7163 

.1354 

.9435 x 10 -1 


0 

.3725 X lO' 3 

.1172 

.8794 X lO -2 

.6953 

.9887 X 10' 1 

.7947 x lO - 1 



.3724 X lO' 3 

.1172 

.8790 X lO* 2 

.6949 

.9884 X 10" 1 

.7944 x lO -2 


-2 

.3392 X 10" 2 

.3148 

.78 57 X 10“ 2 

.6443 

.7967 X lO -2 

.2172 X 10" 1 



.3391 X 10" 2 

.3148 

.7853 X 10“ 2 

.6439 

.7970 X lO -2 

.2172 X 10' 1 


-4 

.3271 X 10" 1 

.3387 

.7608 X lO -2 

.6185 

.9526 X 10 -4 

.2327 X lO* 2 



.3270 X 10" 1 

.3388 

.7604 X lO" 2 

.6182 

.9531 X lO" 4 

.2327 X 10" 2 


-6 

.2576 

.3017 

.6730 X lO -2 

.4338 

.8545 X 10“ 6 

.1846 X 10" 3 



.2575 

.3017 

.6727 X 10~ 2 

.4335 

.8545 X 10" 6 

.1845 X 10 -3 

6000 

1 

0.1189 X 10" 1 

0.2119 

0.8298 X lO -2 

0.6376 

0.3000 X lO' 1 

0.1003 



.1189 X 10" 1 

.2118 

.8294 X lO -2 

.6372 

.2997 X 10 -1 

.1002 


0 

.3597 x 10' 1 

.2937 

.7803 X lO’ 2 

.6122 

.6049 X lO -2 

.4418 X 10' 1 



.3597 X 10" 1 

.2937 

.7798 X lO -2 

.6118 

.6049 X lO -2 

.4416 X 10' 1 


-2 

.2754 

.2955 

.6677 X lO -2 

.4185 

.7144 X 10‘ 4 

.3969 X lO -2 



.2753 

.2951 

.6672 X lO -2 

.4181 

.7131 X 10‘ 4 

.3964 X lO" 2 


-4 

.7352 

.2190 

.4885 X lO -2 

.4078 X 10" 1 

.5370 X 10 -6 

.1075 X 10 -3 



.7338 

.2187 

.4878 X lO -2 

.4062 X HT 1 

.5355 X 10‘ 6 

.1071 X 10" 3 

8000 

2 

0.3655 X 10' 1 

0.1841 

0.8312 X lO -2 

0.6083 

0.2595 X 10 -1 

0.1368 



.36 54 X 10" 1 

.1834 

.8312 X lO -2 

.6079 

.2574 X lO’ 1 

.1362 


0 

.2886 

.2806 

.6684 X 10 -2 

.4056 

.6532 X lO' 3 

.1788 X lO - 1 



.2883 

.2800 

.6677 X 10 -2 

,4048 

.6504 X lO' 3 

.1782 X 10" 1 
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TABLE ID.- COMPOSITION COMPUTED BY PRESENT METHOD AND COMPOSITION FROM REFERENCE 6 - Continued 


[First value of paired numbers, present method; second value, ref. 6] 


(b) Model B 


T, °K 


N 

O 

A 

n 2 

NO 

e" 

N+ 

o + 

4 000 

-3 

(a) 

(a) 

(a) 

(a) 

(a) 

(a) 

(a) 

(a) 



0.1056 X 10' 1 

0.3371 

0.7715 x 10" 2 

0.6360 

0.7372 X lO’ 2 

0.1065 X lO' 4 

0.3972 X lO -10 

0.3768 x 10“ 8 


-4 

(a) 

(a) 

(a) 

(a) 

(a) 

(a) 

(a) 

(a) 



.3270 X 10- 1 

.3388 

.7604 X 10' 2 

.6182 

.2327 X lO' 2 

.1879 X 10“ 4 

.6864 X 10“ 9 

.2115 X 10“ 7 


-6 

(a) 

(a) 

(a) 

(a) 

(a) 

(a) 

(a) 

(a) 



.2575 

.3017 

.6727 X 10~ 2 

.4335 

.1845 X 10“ 3 

.5015 X 10“ 4 

.1792 x 10“ 6 

.6241 X 10" 6 

6000 

-1 

(a) 

(a) 

(a) 

(a) 

(a) 

(a) 

(a) 

(a) 



0.1056 

0.3147 

0.7374 X 10“ 2 

0.5563 

0.1467 X lO -1 

0.1742 X 10" 3 

0.5502 x 10 -6 

0.1984 X 10" 5 


-2 

.2740 

.3012 

.6617 X 10-2 

.4141 

.4026 X lO- 2 

.7200 X 10" 4 

.3091 X 10“ 4 

.4109 X 10 -4 



.2753 

.2951 

.6672 X 10“ 2 

.4181 

.3964 X lO' 2 

.2796 X 10 -3 

.7998 X 10 -5 

.1037 X 10 -4 


-4 

.7342 

.2186 

.4881 X 10-2 

.4066 X 10" 1 

.1071 X lO' 3 

.7679 X lO' 3 

.5645 X lO -3 

.2034 X lO- 3 



.7338 

.2187 

.4878 X 10“ 2 

.4062 X HT 1 

.1071 X lO" 3 

.8662 X lO" 3 

.5002 X lO" 3 

.1804 x lO" 3 


-6 

.7725 

.2072 

.4660 X 10“ 2 

.4717 X 10‘ 3 

.1118 X 10“ 5 

.7591 X lO -2 

.5732 X lO' 2 

.1859 X lO" 2 



.7721 

.2074 

.4650 X 10- 2 

.4712 X 10“ 3 

.1119 x 10“ 5 

.7649 X lO -2 

.5685 X lO" 2 

.1847 X lO" 2 

8 000 

0 

(b) 

(b) 

(b) 

(b) 

(b) 

(b) 

(b) 

(b) 



0.2883 

0.2800 

0.6677 X 10“ 2 

0.4048 

0.1782 x KT 1 

0.7408 X 10" 3 

0.5744 X 10“ 4 

0.4439 X 10" 4 


-2 

.7358 

.2166 

.4852 X lO- 2 

.3624 X 10- 1 

.4835 X 10" 3 

.3021 X lO" 2 

.2448 X lO" 2 

.5732 X lO" 3 



.7354 

.2166 

.4843 X 10" 2 

.3620 X 10" 1 

.4833 x 10“ 3 

.3190 X lO' 2 

.2317 x lO -2 

.5430 X 10" 3 


-4 

.7378 

.1994 

.4560 X 10 -2 

.3878 X lO" 3 

.4750 x 10" 5 

.2891 X lO’ 1 

.2379 X 10" 1 

.5115 X lO" 2 



.7375 

.1996 

.4493 x 10" 2 

.3874 x 10“ 3 

.4753 X 10” 5 

.2899 X 10" 1 

.2371 x 10" 1 

.5105 x 10“ 2 


-6 

.4554 

.1293 

.3728 x lO" 2 

.1808 X 10" 5 

.2327 X 10" 7 

.2058 

.1679 

.3792 X 10“ 1 



.4552 

.1294 

.3202 X lO" 2 

.1806 X 10" 5 

.2328 X lO -7 

.2061 

.1675 

.3790 X 10“ 1 

10 000 

0 

0.6203 

0.2319 

0.5321 X lO" 2 

0.1304 

0.6151 X lO -2 

0.2937 X lO" 2 

0.2364 X 10“ 2 

0.5735 X 10“ 3 



.6198 

.2312 

.5308 X 10“ 2 

.1302 

.6127 X lO -2 

.3470 X lO -2 

.1999 X lO -2 

.4839 X 10“ 3 


-2 

.7373 

.2010 

.4575 X lO" 2 

.2139 

.7360 X 10’ 4 

.2748 X 10 _1 

.2335 X KT 1 

.4130 X 10“ 2 



.7368 

.2011 

.4478 X 10“ 2 

.2136 

.7359 X 10" 4 

.2761 X KT 1 

.2323 X 10“ 1 

.4113 X 10“ 2 


-4 

.4656 

.1376 

.3771 X lO" 2 

.1035 x 10“ 4 

.3859 X 10" 6 

.1965 

.1649 

.3160 x KT 1 



.4654 

.1377 

.3043 x 10“ 2 

.1034 x 10“ 4 

.3861 X 10- 6 

.1969 

.1645 

.3157 X 10 -1 


-5 

.1996 

.6646 X 10“ 1 

.2978 x lO -2 

.2408 x 10' 6 

.1012 x 10" 7 

.3655 

.2984 

.6711 X lO" 1 



.1995 

.6788 X 10" 1 

.1483 x lO" 2 

.2406 x 10 -6 

.1033 x lO -7 

.3654 

.2982 

.6582 X 10“ 1 

12 000 

0 

0.7276 

0.2099 

0.4770 X 10“ 2 

0.2809 x KT 1 

0.2077 x 10’ 2 

0.1382 X 10" 1 

0.1180 X 10" 1 

0.2021 x lO- 2 



.7268 

.2095 

.4703 X lO" 2 

.2803 X 10" 1 

.2076 x lO -2 

.1423 X 10 _1 

.1146 X lO -1 

.1959 x 10 “ 2 


-2 

.6060 

.1720 

.4184 x lO" 2 

.2219 X 10 -3 

.1619 X 10 -4 

.1088 

.9313 X 10" 1 

.1567 X 10 -1 



.6057 

.1721 

.3674 X lO" 2 

.2217 X 10“ 3 

.1619 X 10 -4 

.1091 

.9283 X 10 -1 

.1564 X 10' 1 


-3 

.3807 

.1194 

.3529 X lO” 2 

.1039 X KT 4 

.8368 X 10" 6 

.2482 

.2093 

.3890 X KT 1 



.3805 

.1195 

.2359 X lO" 2 

.1038 X 10" 4 

.8373 X 10“ 6 

.2488 

.2087 

.3885 X 10" 1 

14 000 

0 

0.7088 

0.1970 

0.4537 X lO -2 

0.6840 X lO -2 

.8196 X KT 3 

0.4102 X 10" 1 

0.3535 X 10' 1 

U.5674 X 10“ 2 



.7080 

.1967 

.4305 X 10“ 2 

.6825 X 10 -2 

.8175 x 10“ 3 

.4133 X 10" 1 

.3505 X 10“ 1 

.5623 X lO” 2 


-2 

.3898 

.1224 

.3557 X lO -2 

.2637 X 10“ 4 

.3570 X 10“ 5 

.2421 

.2050 

.3715 X 10 -1 



.3896 

.1225 

.2263 X 10" 2 

.2635 X 10 -4 

.3571 X 10 -5 

.2428 

.2043 

.3707 X 10“ 1 


a Wrong root (see text). 
b No root. 


43 


TABLE EL- COMPOSITION COMPUTED BY PRESENT METHOD AND COMPOSITION FROM REFERENCE 6 - Concluded 
[First value of paired numbers, present method; second value, ref. 6] 

(c) Model C 


r,°K 

io “» (S) 

N 

O 

A 

e~ 

N+ 

o + 

A + 

6000 

-6 

0.7731 

0.2071 

0.4650 X 10" 2 

0.7596 x lO* 2 

0.5732 x lO' 2 

0.1857 x 10" 2 

0.7997 x 10" 5 



.7721 

.2074 

.4650 x 10-2 

.7649 x lO' 2 

.5685 X lO -2 

.1847 x 10' 2 

.7942 x 10" 5 

8000 

-4 

0.7383 

0.1993 

0.4494 x 10” 2 

0.2894 x 10“ 1 

0.2377 x lO' 1 

0.5106 x lO' 2 

0.6409 x 10‘ 4 



.7375 

.1996 

.4493 x 10' 2 

.2899 x 10' 1 

.2371 x 10 -1 

.5105 x lO' 2 

.6398 x 10‘ 4 


-6 

.4554 

.1293 

.3204 X 10‘ 2 

.2060 

.1676 

.3787 X 10* 1 

.5226 x 10" 1 



.4552 

.1294 

.3202 x lO- 2 

.2061 

.1675 

.3790 x lO” 1 

.5221 X lO -3 

10000 

-3 

0.6558 

0.1818 

0.4044 x lO" 2 

0.7916 x 10* 1 

0.6686 x 10" 1 

0.1203 x 10- 1 

0.2777 X lO -3 



.6553 

.1820 

.4043 x lO* 2 

.7920 x 10“ 1 

.6677 x lO” 1 

.1203 x 10* 1 

.2775 X 10“ 3 


-4 

.4657 

.1375 

.3045 x lO -2 

.1969 

.1646 

.3154 x 10" 1 

.7249 X 10" 3 



.4654 

.1377 

.3043 x 10-2 

.1969 

.1645 

.3157 x 10- 1 

.7244 x lO' 3 


-6 

.3913 x 10- 1 

.1539 x 10* 2 

.3323 x lO' 3 

.4726 

.3748 

.9564 x 10* 1 

.2143 x 10' 2 



.3911 x 10 _1 

.1541 x 10' 2 

.3321 x 10' 3 

.4726 

.3746 

.9574 x 10" 1 

.2142 x 10' 2 

12000 

-2 

0.6063 

0.1719 

0.3675 x 10“ 2 

0.1091 

0.9294 x 10" 1 

0.1563 x 10' 1 

0.5068 x 10‘ 3 



.6057 

.1721 

.3674 x 10“ 2 

.1091 

.9283 x 10" 1 

.1564 X 10" 1 

.5066 x lO -3 


-4 

.1299 

.4927 x 10" 1 

.8370 x 10" 3 

.4100 

.3331 

.7493 x 10" 1 

.1932 x 10' 2 



.1298 

.4932 x 10" 1 

.8365 x 10" 3 

.4100 

.3330 

.7501 x 10‘ 1 

.1931 x 10' 2 


-6 

.2212 x 10" 2 

.9970 x 10" 3 

.1470 x 10" 4 

.4984 

.3914 

.1046 

.2340 x 10~ 2 



.2211 x 10" 2 

.9980 x 10" 3 

.1468 x 10“ 4 

.4984 

.3912 

.1047 

.2337 x 10‘ 2 

14000 

-1 

0.6084 

0.1728 

0.3597 x 10" 2 

0.1076 

0.9195 x 10" 1 

0.1508 x 10" 1 

0.5912 x 10" 3 



.6074 

.1729 

.3598 x 10" 2 

.1076 

.9181 x 10“ 1 

.1509 x 10" 1 

.5914 x 10‘ 3 


-2 

.3898 

.1224 

.2264 x 10" 2 

.2428 

.2044 

.3704 x 10* 1 

.1291 x 10' 2 



.3896 

.1225 

.2263 x 10” 2 

.2428 

.2043 

.3707 x 10- 1 

.1290 x 10‘ 2 


-4 

.2287 x 10“ 1 

.1021 x 10“ 1 

.1264 x 10" 3 

.4834 

.3826 

.9855 x 10" 1 

.2298 x 10‘ 2 



.2286 x 10" 1 

.1022 x 10" 1 

.1263 x 10“ 3 

.4834 

.3823 

.9865 x 10“ 1 

.2296 X 10‘ 2 


-6 

.2539 x 10" 3 

.1180 x 10" 3 

.1397 x 10" 5 

.4998 

.3923 

.1052 

.2346 X 10‘ 2 



.2533 x 10“ 3 

.1181 x 10" 3 

.1357 x 10" 5 

.5002 

.3910 

.1052 

.2278 x 10' 2 


(d) Model D 


r, °K 

log i°(£) 

e~ 

N + 

o+ 

A + 

N ++ 

o++ 

A++ 

13 000 

-6 

0.5001 

.4996 

0.3922 

.3917 

0.1052 

.1051 

0.2336 X lO -2 
.2333 x lO' 2 

0.1060 x 10' 3 
.1059 x 10' 3 

0.6773 x 10' 6 
.6770 x lO -6 

0.1027 x lO' 4 
.1026 x 10‘ 4 

14 000 

-6 

0.5004 

.5002 

0.3913 

.3910 

0.1052 

.1052 

0.2279 x 10" 2 
.2278 x 10“ 2 

0.7689 x lO' 3 
.7687 x lO' 3 

0.6957 x 10“ 5 
.6962 x 10" 5 

0.6555 x 10' 4 
.6 554 x 10' 4 

15000 

-5 

0.5002 

.4995 

0.3918 

.3911 

0.1052 

.1050 

0.2311 x 10" 2 
.2308 x 10" 2 

0.4370 x lO' 3 
.4369 X lO’ 3 

0.5314 x 10' 5 
.5312 x 10 -5 

0.3439 x 10' 4 
.3439 x 10‘ 4 


-6 

.5023 

.5022 

.3863 

.3861 

.1047 

.1048 

.2038 x 10" 2 
.2036 x 10" 2 

.4242 x 10' 2 
.4240 x 10" 2 

.5208 x 10‘ 4 
.5213 x 10‘ 4 

.2985 x 10' 3 
.2983 x 10' 3 
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TABLE IV.- 3 d n / 8 rj FOR MODEL A 
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" b 2 8 ( r 0 ' b 2) 

-^[ b 2 ( b 2 +b 3 ) + (> 


2 + b 0 + b^l 

+ h 3 J 

X jl + b2 + bg + 2 Tg ( 

rfbo'T’o + 

4 y n (y n - bo') 

- f*2v-< (bo + b 



TABLE V.- 3d n /arj FOR MODEL C 





TABLE VI.- 9d n /9 T j FOR MODEL D 



NASA.- Langley , 1967 24 L-5544 


"The aeronautical and space activities of the United States shall he 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof ” 

— National Aeronautics and Space Act of 1958 
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tion because of preliminary data, security classification, or other reasons. 
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under a NASA contract or grant and considered an important contribution to 
existing knowledge. 
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SPECIAL PUBLICATIONS: Information derived from or of value to NASA 
activities. Publications include conference proceedings, monographs, data 
compilations, handbooks, sourcebooks, and special bibliographies. 

TECHNOLOGY UTILIZATION PUBLICATIONS: Information on tech- 
nology used by NASA that may be of particular interest in commercial and other 
non-aerospace applications. Publications include Tech Briefs, Technology 
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